Interaction of quasilocal harmonic modes and boson peak in glasses 
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The direct proportionality relation between the boson peak maximum in glasses, u>b, and the Ioffe- 
Regel crossover frequency for phonons, Ud, is established. For several investigated materials ujb = 
(1.5 ± Q.l)u)d- At the frequency uid the mean free path of the phonons I becomes equal to their 
wavelength because of strong resonant scattering on quasilocal harmonic oscillators. Above this 
frequency phonons cease to exist. We prove that the established correlation between ujb and uid 
holds in the general case and is a direct consequence of bilinear coupling of quasilocal oscillators 
with the strain field. 
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The properties of harmonic vibrational excitations in 
disordered media and glasses become now a very active 
topic of scientific research . Contrary to the quite well 
established behavior of electrons in disordered conduc- 
tors there is no consensus at all regarding the harmonic 
vibrations in ordinary glasses. The most common and 
challenging of their signatures is the so-called boson peak 
observed in numerous experiments in the low-frequency 
Raman and inelastic neutron scattering. The physical 
origin of the peak is however still a matter of great de- 
bates. The common view is that the solution of this prob- 
lem is a corner-stone for our veritable understanding of 
glassy vibrational dynamics. 

The main discussion in the literature involves now the 
question whether the harmonic vibrations responsible for 
the boson peak are propagating plane waves (phonon 
like) H or localized because of disorder jjj. The third 
possibility, which we share in this paper, is that they 
are neither propagating waves nor localized but have a 
diffusive nature Q]. 

For the solution of this crucial question quite power- 
ful and expensive experimental techniques are now in use. 
First of all these are Raman experiments themselves. But 
since visible light due to energy and momentum conser- 
vation laws does not interact properly with sound-like ex- 
citations, one has to use two other possibilities, namely 
inelastic X-ray and neutron scattering. 

The main difficulty in the inelastic X-ray scattering ex- 
periments is the very high incident photon energy (~ 20 
KeV), and its relatively small change, of the order of 
1CP 7 . As a result, the "X-ray boson peak" is superim- 
posed on the steep wings of strong elastic line in the 
forward direction. Therefore, a serious problem arises 
to correctly resolve it and answer the crucial question 
whether it changes with the momentum transfer. 

This difficulty explains why recently two independent 
groups using the same experimental setup but different 
fitting procedures have arrived to completely opposite 
conclusions about the propagating character of the exci- 



tations in vitreous SiC>2 at the boson peak range |^,^|. 
The situation with inelastic neutron Brillouin scattering 
(for small momentum transfer) is not much better. 

To solve this problem from our point of view it is neces- 
sary to separate the Brillouin (phonon) lines of the boson 
peak (as in the usual light scattering experiments). Then 
increasing the momentum transfer q one can observe how 
the Brillouin line shifts and broadens approaching the 
boson peak from the low frequency side. This scenario 
implies the existence in the q range 0.1 -j- 1 nm _1 of two 
peaks, one from damped phonons (Brillouin line) and an- 
other the boson peak itself. However such experiments 
will hardly appear in a nearest future. 

Therefore at this stage of our knowledges any theoreti- 
cal insight into the problem would help a lot to establish 
the true picture of the harmonic vibrations in glasses. 
The main idea is that the boson peak, being a universal 
glassy property, should be related with other universal 
properties of glasses. 

In the papers J5|-|7j the idea was put forward that the 
boson peak corresponds to the Ioffe-Regel crossover fre- 
quency for phonons, L04. As it was conjectured in Q] 
above ujd phonons cease to exist as well defined plane 
wave excitations. Later this idea was adopted (sup- 
ported) in many papers ||,|| , though taken alone without 
any theory, it of course could not explain why actually 
this correlation occurs. 

The explanation of such a correlation based on the 
theory of soft atomic potentials in glasses § was pro- 
posed in |t|. The main idea is that at this frequency the 
dipole-dipole interaction between quasilocal harmonic os- 
cillators renormalizes their density of states (DOS) from 
the bare value g(u>) oc w 4 (for independent oscillators) 
to the g(uj) oc to (for coupled oscillators) behavior. As a 
result a boson peak appears in g(uj)/uj 2 at a frequency 
~ ujd- The reason for the DOS transformation is a strong 
level repulsion in the course of the interaction of oscilla- 
tors. Above uid the oscillators become delocalized pfj| . 

This correlation was checked indirectly for As x Sei— x 
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glasses || and excellent agreement was found between 
the position of the bump in the reduced specific heat, 
C(T)/T 3 (prototype of the boson peak), and the energy 
E c i = hud for different x compositions. The direct evi- 
dence should of course include the comparison of the bo- 
son peak and the Ioffe-Regel crossover frequency which is 
explicitly calculated from the mean free path of phonons 
I. In its turn I can be easily evaluated from the fit of the 
thermal conductivity data which are available for many 
glasses. The purpose of our paper is just to make such a 
comparison. 

In this paper on the basis of experimental data analysis 
for different glasses we have established for a first time 
a direct proportionality relation between the position of 
the boson peak and the Ioffe-Regel crossover frequency 
LJd- At this frequency the phonon mean free path with 
respect to the resonant scattering on quasilocal harmonic 
oscillators becomes equal to the wave length. To proceed 
further we should briefly remind the main theoretical and 
experimental results in this field. 

It is well known that for many low-temperature prop- 
erties of glasses two-level systems (TLS's) and phonons 
are responsible However, at higher frequencies it 

was well established that there is an excess of additional 
low frequency harmonic modes which dominate the spe- 
cific heat above a few Kelvin [ p~2|JT^ ] . They also strongly 
reduce the phonon mean free path producing the plateau 
in the thermal conductivity of glasses [jl4j . 

In Ref. § it was argued that these excess modes are 
quasilocal soft harmonic oscillators (HO) coexisting with 
TLS's (and phonons) with density of states increasing 
as g(ui) oc uo A at moderate frequencies {W/h < ui < uod, 
see below). As a result with rising the temperature the 
linear temperature behavior of the specific heat (TLS's 
contribution) changes to C(T) oc T 5 dependence. Such 
crossover results in a minimum in C(T)/T 3 , at some 
T min «0.5-3K |p. 

The quasilocal vibrations open also a new effective 
channel for the phonon scattering. Their steeply rising 
DOS, g(uj), leads to the same frequency dependence of 
the inverse mean free path for phonons because of reso- 
nant scattering on these oscillators p3 
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This dependence looks like elastic Rayleigh scattering 
from glass inhomogeneities though the physical mecha- 
nism is drastically different. We believe that in ordinary 
glasses the Rayleigh scattering is small in comparison 
with this contribution. 

Except for the sound velocity v there are only two pa- 
rameters in ([[]). The first one is a characteristic energy 
W which is related with the position of the minimum in 
C(T)/T 3 , W « 2fcT min [HJ]. The second is a coupling 
constant C which describes the relative change of v with 
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temperature due to resonance scattering of phonons on 

tls's 



(Av/v) xeStThS = Chx(T/T ) 
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Both experimental values of C and T m j n are well known 
for many glasses. It gave the unique possibility, with- 
out any fitting parameters using ([!]), to successfully re- 
produce the correct value of the thermal conductivity 
plateau for vitreous SW2 and Se jl4| and for some other 
glasses as well [|l5| . From this point of view one can con- 
sider Eq. (0) to be in a good agreement with the thermal 
conductivity data. 

Using ([!]) one can easily obtain the_expression for the 
Ioffe-Regel crossover frequency Ud 
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For example for v - Si0 2 , Q = 3.1 x 1CT 4 , W/k = 4.2 K 
and hujd/k w 47 K. 

Now we can proceed to the main aim of our paper and 
check whether there exists a correlation between the bo- 
son peak frequency uif, and u>d determined by (J3J) . For 
that purpose we collected in the Table | all the neces- 
sary experimental parameters r m i n , Ci and tub for several 
glasses. The two different values of T m i n for a — Ge02 
correspond to two different samples investigated in p8| . 
The value of T m ; n for LaSF-7 glass is a reasonable guess 
value (we did not find it in a literature). 

Now using the relation W = 2kT min and Eq. || we 
can calculate u>d for longitudinal phonons and compare 
it with cub- The result of this comparison is shown on 
Fig.l. One can see from the figure that all the data lies 
near a straight line which has the slope equal to 1.53. It 
means that within the experimental accuracy there is a 
direct proportionality relation between the boson peak 
maximum cob and the Ioffe-Regel crossover frequency for 
phonons LUd- 
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In what follows we are going to prove that this corre- 
lation results in fact from the bilinear coupling between 
quasilocal oscillators and the deformation field e 
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proposed in |14|. In this sense it can be regarded as a 
rather general rule. 

With the use of (|j) one can easily show that the elastic 
interaction between two oscillators has a dipole-dipole 
character 
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with p being the mass density of the glass. 

Suppose that the DOS of noninteracting oscillators is 
g(u>). Then making use of (Q) and the Golden rule one 
can calculate in the usual way the phonon mean free path 
due to the resonance scattering on these oscillators 
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where M is an oscillator mass. It is important that this 
scattering has the same frequency dependence as the bare 
density g(u>). The Ioffe-Regel crossover frequency ujd can 
be estimated from the equation 
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In order to determine a position of the boson peak we 
should consider the interaction between resonant oscilla- 
tors. The concentration of such oscillators with frequen- 
cies u> lying in a small interval 5u> around u> is given by 
the product g(uj)Suj. Therefore the non-diagonal tran- 
sition matrix element of the coupling between the two 
oscillators (||) is equal 



A J = (n i+ i,nj\Vi nt \ni,rij + i) 
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The physical meaning of this quantity is that it gives a 
characteristic value for a level repulsion of two resonant 
oscillators because of their interaction. 

If AJ < USu) the level repulsion is weak and the in- 
teraction does not change g(u>). If A J > USco the level 
repulsion is strong and renormalizes the original density 
0. The new DOS can be found dividing the oscillator 
concentration in the interval by the characteristic value 
of their level repulsion 
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We see that the result is a linear function of lj which 
depends only on material parameters of the glass but 
does not depend on the bare density g(w). 



In the case when the reduced bare DOS g(uj)/ui 2 in- 
creases with u> we will have a boson peak at some 
crossover frequency from weak to strong coupling. This 
frequency is determined by the equation A J = hdu. One 
can easily see that this equation coincides with equation 
(jjj) for ujd- This explains why ujd correlates with the po- 
sition of the boson peak in glasses. 

The resulting linear u> dependence of vibrational DOS 
above the boson peak is a very important and general fea- 
ture of a system of interacting harmonic oscillators with 
bilinear dipole-dipole coupling (||). It is a direct conse- 
quence of linear u> dependence of | a^TT, i | 2 = 2Mu>/h 
and therefore can be regarded as universal. It holds for 
any bare g{u>) of quasilocal harmonic modes which rises 
with frequency faster then to 2 . It was observed in ex- 
periments on inelastic neutron scattering for amorphous 
polymers ||, a - Si0 2 @ and a - GeSe 2 @. 

The picture described above implies derealization of 
the harmonic excitations above the boson peak though 
they are not supposed to be plane waves. This is also 
consistent with existing numerical results. The fact that 
the Ioffe-Regel criterion does not mean localization for 
vibrational excitations in disordered system was acknowl- 
edged in the computer analysis of the percolation network 
p5| . It was found that vibrational modes above the Ioffe- 
Regel crossover give substantial contribution to the heat 
transfer and have some diffusive nature. 

Similar results were obtained in [ p6[ for a — S1O2 where 
all modes with Tito between 5 — 110 meV were found to 
be delocalized (just to the right of the boson peak). For 
amorphous silicon it was found that only 3% of vi- 
brational modes (above 70 meV) are localized and do 
not participate in the heat transfer. The remaining 97% 
of the modes are delocalized and the majority of them 
(93%) are not plane waves. 

In numerical experiments |3q| the truly localized modes 
were discovered at the very end of the frequency spec- 
trum. It is worthwhile to mention that at low frequencies 
the calculated DOS has a part with linear lj dependence 
just in the region of delocalized states. The same results 
were obtained for amorphous selenium p9| . The modes 
with frequencies above 0.5 THz were delocalized (again 
just above the boson peak for this glass). 

To conclude, we have established for the first time 
the correlation between the boson peak frequency and 
the Ioffe-Regel crossover frequency. We proved that this 
correlation occurs because of the bilinear coupling be- 
tween sound waves and quasilocal harmonic vibrations 
in glasses. 
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Figure caption 
Fig.l Position of the boson peak ojb versus the Ioffc- 
Regel crossover frequency Ud for several glasses. 
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